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The deconfinement transition in 5(7(4) lattice gauge theory is studied on Nf x Nt lattices with Ns = 8-16 and 
Nt = 4-8 using a modified Wilson action which is expected to have no bulk transitions. The susceptibility x^{l\^ 
is found to increase linearly with spatial volume for Nt — 4, 5, and 6, indicating a first order deconfinement phase 
transition. The latent heat is estimated to be « | of the corresponding ideal gas energy density at Tc. 



1. INTRODUCTION 

The nature of the phase transition(s) to Quark- 
Gluon Plasma, which we hope to see in the ex- 
periments of RHIC and later of LHC, and the 
physics driving it, have always been of great in- 
terest. While the real world has presumably two 
very light {u,d) flavours of quarks and one some- 
what heavier (s) flavour, both analytical and nu- 
merical methods in lattice QCD begin from the 
limiting cases of either massless or infinitely mas- 
sive quarks. One talks of the chiral symme- 
try restoring phase transition and the deconfine- 
ment phase transition in these two cases respec- 
tively and has suitable order parameters to in- 
vestigate them. For quarks with N colours and 
Nf massless flavours, these transitions are related 
to spontaneous breaking of a global Z{N) and 
SU{Nf ) X SU{Nf) chiral symmetry. Since these 
symmetries are broken explicitly to various ex- 
tents in the real world, which of them is more 
relevant is a priori not clear. The low masses of 
the light flavours suggest chiral symmetry to be 
the dominant one. However, it is seen in numer- 
ical simulations that the energy density shows a 
large change at the chiral transition and even the 
order parameter for the deconfinement transition. 



1 ^' 
t=i 



'E-mail: gavai@tifr.res.in 



(1) 



also rises to nonzero values there. These appar- 
ently mysterious observations can be explained [0 
using large N arguments, if the deconfinement 
transition for > 4 is of second order. SU (4) is 
clearly the simplest case to check this out. 

Numerical simulations of SU (4) theory at finite 
temperatures have been done in the pastj2j and 
recently 1^ as well. All of them used the Wilson 
action, or the more general mixed action: 
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A well known problem in the simulations with 
these actions, especially for large N, is the pres- 
ence of a bulk transition which is a lattice arti- 
fact. The phase diagram of the mixed action. 
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Figure 1. A schematic phase diagram in {P,Pa)- 
plane for the mixed action of eq. (||). 

in its coupling plane is as shown in Fig. ^ The 
solid lines in it show first order bulk transitions 
lines. The dotted line after the end point D is 
drawn to suggest the impact D may have on the 
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Wilson axis (/3a = 0). For > 4, D is expected 
to be where E is shown, causing a first order bulk 
transition for the usual Wilson action. In order 
to avoid it, simulations at negative Pa H] and for 
larger ||^ A^f = 6 were made. They lead to a first 
order deconfinement phase transition for SU (4) . 

From our extensive studies Q of the deconfine- 
ment phase transition for the action above but 
for the SU{2) theory, we know that bulk tran- 
sitions affect the order and location of the de- 
confinement transition in subtle and inexplicable 
ways, even leading to apparent qualitative viola- 
tions of universality. Universality was restored ||^ 
in that case only after eliminating the bulk tran- 
sitions associated with the Z{2) vortices and Z{2) 
monopoles by adding large chemical potentials for 
them. It seems natural to expect that the bulk 
transitions for N > 2 can also be cleaned off by 
suppressing the corresponding Z{N) objects. We 
pursue this idea here for 5't/(4) to investigate its 
deconfinement phase transition. 

2. SIMULATIONS AND RESULTS 

Generalizing the idea of positive plaquette 
modelsQ in the literature, we use the action 
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where — tt < a < tt is the phase of tr Up. By 
adding the adjoint term of eq. (2) to the action 
(3), one sees that the phase diagram of the resul- 
tant mixed should not have the bulk lines AC or 
BC and hence the endpoint D or E. 

We have simulated the above action on x Nt 
lattices for Ns =8, 10, 12, 15, 16 and Nt = 4, 
5, 6, 8 using a 15-hit Metropolis et al. algo- 
rithm. The calculations were done on a cluster 
of Pentiums. Typically short runs to look for 
points of rapid variations in were followed by 
long runs (a few million sweeps) to determine the 
susceptibility x\l\ using the histogramming tech- 
nique. Usual finite size scaling techniques were 
used to determine the order of the transition and 
its exponents. 

In simulations on x 4 lattices, Ng — 8, 10, 
12, one sees hot and cold starts converge quickly 
at couplings a little away from the transition 



point on its both sides but a clear co-existence 
of states is visible for all lattices at the transition 
point. The tunneling frequency goes down with 
spatial volume. The histograms of \L\ show peaks 
which become narrower with increasing volume 
and the gap between them remains unchanged. 
These classic signs of a first order phase transi- 
tion are confirmed by a quantitative analysis of 
the linear growth of x\l\ with volume, as seen in 
Fig. H The horizontal lines in each case are pre- 
dictions obtained by scaling the A'^, = 8 results 
linearly with volume. 
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Figure 2. The susceptibility x\l\ as a function 
of /? for X 4 lattices. 

For larger Nt , we used many longer runs in the 
region of strong variation of to obtain the 
susceptibility directly and used the histogram- 
ming only for the finer determination of the criti- 
cal coupling. Our results for (|L|) as a function of 
(3 clearly show the expected shift for a deconfine- 
ment phase transition for Nt — 5, 6 and 8. This is 
evident in the f3c determinations from the x|l|'^i 
as seen in Table 1 for Nt = 6 for two different 
spatial volumes. Again, using the peak height 
for the smaller spatial volume, the x"^'^^ on the 
bigger lattice can be predicted, assuming a lin- 
ear growth with volume. The prediction listed in 
the table can be seen to be in very good agree- 
ment with the direct Monte Carlo determination. 
Along with the shifts in f3, it confirms that the 
same physical deconfinement phase transition is 
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Table 1 

The values of (3 at which long simulations were 
performed on x 6 lattices, /3c and the height 
of the |L|-susceptibility peak, xTl\^- 





/3 






^max 
Apredcted 


12 


10.675 


10.686(5) 


4.36(35) 




16 


10.675 


10.676(5) 


10.43(95) 


10.3(8) 



being simulated on these lattices thus approach- 
ing the continuum hmit of a ^ in a progressive 
manner by keeping the transition temperature Tc 
constant in physical units. 



3. LATENT HEAT 

While the results above for the deconfinement 
order parameter and the corresponding sus- 
ceptibility, are indicative of a first order de- 
confinement phase transition, one needs to make 
sure that they indeed are not due to a coincident 
first order bulk transition. Apart from the charac- 
teristic (logarithmic) shift of the transition point 
with iVt, seen above, the latent heat of a first 
order deconfinement phase transition should also 
remain constant as Nt oo. Requiring the pres- 
sure to be continuous at the deconfinement phase 
transition, the latent heat can be obtained from 
two different observables Ai = A{e — 3p)/T^, and 
A2 = A{e+p)/T^, where 

oa 

A2 = 32iV4 (APt-APs) , (4) 

andC(g2) = (1-0.236652-^0(54)) for S';7(4). A 
denotes discontinuities across the transition in the 
respective variables. In order to obtain the AP, 
APs and APt, the minimum of the histogram 
N{\L\) was used to separate the two phases in 
each case. From eq. (jj), it is clear that the 
plaquette discontinuity AP oc N^* in order to 
obtain the same latent heat in physical units, as 
Nt 00. Indeed, its decrease with Nt was seen 
to be consistent with expectations for Nt > 5. 
Furthermore, both estimates must agree in this 
limit, as the neglected cut-off corrections become 
then insignificant. Table suggests this to be the 



Table 2 

Both the latent heat estimates of eq.(||) as a func- 
tion of Nt- 



Nt 


4 


5 6 


8 


Ai 


21.03(5) 


11.02(6) 8.31(5) 


6.57(16) 


A2 


9.89(14) 


7.77(40) 6.04(60) 


6.45(99) 



case, leading to an estimate which is w | of the 
ideal gas energy density at Tc and agrees with 
earlier results ||]. 



4. SUMMARY 

The linear growth of xfE]^ with volume for Nt 
= 4, suggests a first order deconfinement phase 
transition for SU{A). Various indicators, such 
as, histograms and evolutions, are in accord with 
this. Increasing Nt to 5, 6 and 8 shows the ex- 
pected shift of the deconfinement transition, with 
a growth in x'\l\^ that is consistent with being lin- 
ear in volume. The plaquette discontinuity AP 
decreases as the fourth power of iVj, indicating 
both a lack of a bulk transition and a first order 
deconfinement phase transition. The large esti- 
mated latent heat, being about 2/3 of the cor- 
responding ideal gas energy density, suggests the 
deconfinement transition to grow stronger in na- 
ture as the number of colours in increased. 
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